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Abstract 

We investigate a U(l) gauge invariant quantum mechanical system on a 2D noncom- 
mutative space with coordinates generating a generalized deformed oscillator algebra. 
The Hamiltonian is taken as a quadratic form in gauge covariant derivatives obeying 
the nonlinear Dolan-Grady relations. This restricts the structure function of the de- 
formed oscillator algebra to a quadratic polynomial. The cases when the coordinates 
form the su(2) and sl(2,R) algebras are investigated in detail. Reducing the Hamil- 
tonian to ID finite-difference quasi-exactly solvable operators, we demonstrate partial 
algebraization of the spectrum of the corresponding systems on the fuzzy sphere and 
noncommutative hyperbolic plane. A completely covariant method based on the notion 
of intrinsic algebra is proposed to deal with the spectral problem of such systems. 
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1 Introduction 



The Dolan-Grady relations (DGR) arose 20 years ago in a very interesting paper pp, where 
Dolan and Grady considered a class of systems described by a Hamiltonian given as a linear 
combination of the two operators obeying some nonlinear relations. The relations, afterward 
called by their names, guarantee the existence of an infinite set of mutually commuting 
charges which includes the Hamiltonian. This naturally connects the construction to the 
integrable systems [TJ 12 El HI E] • A rich algebraic structure associated with the Dolan-Grady 
relations was revealed by Perk [2] and Davies [3] who showed that the two operators satisfying 
DGR generate recursively the infinite-dimensional Onsager algebra. The latter appeared in 
a seminal work jH] on exact solution of the planar Ising model. The set of commuting charges 
constructed by Dolan and Grady is linear in generators of the Onsager algebra. 

Recently it was observed |H1 E] that the Dolan-Grady relations arise as necessary con- 
ditions for the anomaly-free quantization of pseudo-classical systems with nonlinear holo- 
morphic supersymmetry (n-HSUSY). In accordance with general algebraic construction of 
the n-HSUSY jHj, the Hamiltonian depending on a natural parameter n is a quadratic form 
in two mutually conjugate operators subjected to DGR, whereas the supercharges are (anti) 
holomorphic polynomials of the order n in the same operators. It was shown that the n- 
HSUSY Hamiltonian belongs to an infinite set of mutually commuting even operators being 
quadratic in the Onsager algebra generators, and that these even integrals together with the 
odd supercharges form a polynomial superalgebra of the order n. Realization of this con- 
struction in various physical systems [E1EI1E] revealed an intimate relation of the n-HSUSY 
with the quasi-exactly solvable (QES) systems [T2] . 

Though the quadratic set of commuting integrals admits any real and even complex values 
for the parameter n playing the role of a coupling constant, such an extension breaks the 
nonlinear supersymmetry. Nevertheless, in Ref. jTJ it was observed that the quasi exact 
solvability of the systems associated with the quadratic set of integrals of motion roots not in 
the n-HSUSY construction itself, but rather in the Dolan-Grady relations. In particular, it 
was shown that the corresponding systems can be quasi-exactly solvable not only for integer 
values of the coupling constant n. 

In the present paper, exploiting the Dolan-Grady relations we consider a class of U(l) 
gauge invariant quantum mechanical systems on a two-dimensional noncommutative mani- 
fold, and demonstrate the quasi exact solvability of their spectral problem for any value of 
the coupling constant. 

The layout of the paper is as follows. In Section 2 the basic information on the Dolan- 
Grady relations and Onsager algebra necessary for further analysis is presented in the con- 
text of integrable systems, n-HSUSY and low-dimensional matrix models. Realization of 
the Dolan-Grady relations in terms of the generalized deformed oscillator algebra is given 
in Section 3. In Section 4 we construct a U(l) gauge invariant quantum mechanical system 
with noncommutative coordinates generating the generalized deformed oscillator algebra. 
In Section 5 we demonstrate that by a reduction to finite-difference equations, the spectral 
problem of the model on the fuzzy sphere and noncommutative hyperbolic plane can be par- 
tially solved for some values of the coupling constant. In Section 6 we develop an algebraic 
covariant approach to deal with the spectral problem, and show that really our noncommu- 
tative system is quasi-exactly solvable for any value of the coupling constant. The obtained 
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results are summarized in Section 7. In Appendix we discuss a deformation of the s[(2,R) 
scheme associated with the finite-difference quasi-exactly solvable equations. 

2 Dolan-Grady relations 

2.1 Dolan-Grady relations and Onsager algebra 

The Dolan-Grady relations pQ can be represented in the form [TOl 

[Z, [Z, [Z, Z]}}=cu 2 [Z, Z] , [Z, [Z, [Z, Z]}}=u 2 [Z, 2] . (2.1) 

Here Z and Z are some operators which we shall call generating elements while oj and u are 
some constants. The nonlinear relations (|2.1jl can be used to construct the infinite set of 
mutually commuting operators for some integrable systems. In the most elegant way they 
can be represented in terms of the Onsager algebra [6 spanned by the generators A n and 
G n having the commutation relations 

[A m , A n ] = 4:G m - n , [G m , A n ] = 2A n j r , m — 2A n _ m , [G m , G n ] = 0, (2-2) 

where m, n e Z. In fact, the operators Aq ~ Z and A\ ~ Z are the generating elements 
of the Onsager algebra obeying the Dolan-Grady relations and all the other operators are 
defined recursively by the relations (|2.2jl |3J. The values of the constants u and Q can be 
fixed by rescaling the generators. Then, the infinite set of the mutually commuting operators 
constructed by Dolan and Grady pQ is represented as 

2J m = A m + A_ m + A (v4 m+ i + Ai_ m ) , (2-3) 

where A G M. The commutativity of the operators J m follows directly from the algebra ([2.2)1 . 
Now, treating the operator J = A + XAi as a Hamiltonian, one obtains, generally, an 
integrable system with infinite number of conserved commuting charges. The charges of the 
form ()2.3)1 appear in various models, such as the Baxter eight- vertex jH] , the two-dimensional 
Ising [Oj and Z^ spin models The power of the algebraic formulation (|2.2jl . ()2.3|1 is that 
it does not refer to the number of dimensions or to the nature of the space-time manifold 
which can be lattice, continuum or loop space. 

2.2 Structure of n-HSUSY 

In the systems associated with the linear set of integrals ()2.3|) the operators Z ( ~ A\) and 
Z ( ~ A ) are supposed to be Hermitian. In Ref. 9J it was shown that the set (12.3)1 is not a 
unique set of mutually commuting operators which can be constructed in terms of generators 
of the Onsager algebra. There exists a set of commuting operators which are quadratic in the 
generators of the algebra ()6.2j) . Such commuting operators naturally appear in the scheme 
of the n-HSUSY (HI EH El HI] , which we discuss shortly below. Note, that for the systems 
associated with the quadratic set of commuting operators, the Z, Z and the u, u are usually 
supposed to be mutually conjugate operators and complex constants, respectively. 
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The Hamiltonian of the system associated with the quadratic set of commuting charges 
(see Eq. (16. 2|) for their explicit form) reads as 



^ = \ {Z,Z} - \ [Z, Z] , (2.4) 

where a real parameter A serves as a coupling constant, while mutually conjugate operators 
Z and Z are supposed to obey the nonlinear Dolan-Grady relations (j2.ll) 1 . 

When the coupling constant is integer, A = ±n, n G Z + , the system ()2.4|) can be extended 
to the supersymmetric system with the Hamiltonian TC n = H n a_a + + H_„a + a_, a± = 
|((7i ± i<7 2 ). In this case the Dolan-Grady relations (|2.1j) guarantee the existence of the odd 
integrals of motion (supercharges) 

n— 1 n—1 _ 

Q = a + J[(z-~(2k-n+\j), Q = a J] [z - |(2fc - n + 1)) . (2.5) 

k=0 k=0 

The supercharges generate a nonlinear superalgebra of the order n. Its structure and further 
details of the algebraic construction underlying the nonlinear holomorphic supersymmetry 
can be found in Ref. 0. 

The advantage of the supersymmetric formulation is that the zero modes of the super- 
charge Q (or Q) form an invariant subspace with respect to the action of the Hamiltonian 
H n (or H_ n ). Furthermore, as a rule the zero modes are easier calculable than the eigen- 
functions of the Hamiltonian due to the factorized structure of the supercharges. When the 
number of the zero modes is finite, then the Hamiltonian is quasi-exactly solvable operator. 
The information on the functional form of the zero modes is very helpful for solving the 
corresponding spectral problem. 

Although the system ()2.4j) formally admits an infinite number of integrals, nevertheless, 
for any given realization of the generating elements in a system with finite number of degrees 
of freedom this infinite set is reduced to a finite set of independent integrals. Such systems 
being not integrable in general case, however, reveal a profound connection with quasi-exactly 
solvable systems [§]. We shall return to this aspect below. 

2.3 Matrix models 

Here we show that the Dolan-Grady relations can be revealed in some low dimensional matrix 
models. 

So, let us consider a matrix model given by the action 



S = - Tr Q [X u , X u ] [X", X"} + R^XW^ 



[2.6) 



where X„ are Hermitian operators on a Hilbert space and R^ u are c-numbers. The last term 
can be treated as describing the interaction of the vector field X u with an external symmetric 

principle, one can also deal with Hermitian operators Z and Z. Then the Hamiltonian (and the whole 
set of the integrals l|6. 2(1 ) is Hermitian if the coupling constant is pure imaginary However, here we will 
not discuss such systems. 
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field R^ v . One can also think that its diagonal part contains the mass term of the vector 
field. The matrix models of the form (|2.6|) are widely used in the context of string theories. 
In the case of positively definite metric the action (|2.6|) corresponds to the potential part 
of the BFSS model ^3] and defines static (vacuum) configurations of the theory. On the 
other hand, for the Lorentz metric the action can be interpreted as an effective theory of the 
reduced model of Yang- Mills theory ^3] with the specific interaction. Actually, action ()2.6|) 
represents a slightly modified massive matrix model [To] . 
The action (|2.6|) leads to the equations of motion 

[X", [X„ X u ]]+R,„X» = 0. (2.7) 

In the two-dimensional case Eq. (|2.7|) can be represented in the form 

[X_, X + }} + i? + _X_ + i?__X + = 0, 
[X+, [X+, + R + -X + + i? ++ X_ =0, 

where X± = ^{X\ ± X 2 ) or X± = -^(Xi ± iX 2 ) for the Lorentz or Euclidean metric, 

respectively. In the first case R ++ and R are real constants, while in the latter R*__ = R ++ . 

It is easy to see that any solution of Eq. (|2.8j) also obeys the Dolan-Grady relations (|2.1|) 
with the identification: 

X_ ~ Z, X+ ~ Z, R ~ — cj 2 , R++ ~ — 

Moreover, note that in the case of positively definite metric and when -R++ = R = 

0, the equations ()2.8|) coincide with the basic commutations relations for the single-mode 
parafermion [T5] . 

Now let us consider the following modification of the three-dimensional matrix model 



f\ 4 \ 

S — — Tr I - [Xj, Xk] — -iatjkiXjXkXi + RjkXjXk J 



;2.9) 



where a and Rjk are real numbers and we imply positively definite metric. This model can 
be treated as the matrix model corresponding to the usual Yang-Mills theory with Chern- 
Simon term. As before, the quadratic term in action (j2.9J) represents a specific interaction 
with an external symmetric field. The equations of motion are 

[[X k , Xj], X k ] - iae jkl [X k , X t ] + R jk X k = 0. (2.10) 

Imposing the condition 

[X ,X ± ] = ±PX ± (2.11) 

with (3 G R on the components X , X±, the equation for j = is satisfied identically for 
(3 — a, Roi = 0. The equations of motion for X± have the form ()2.8|) with the substitution 
— > + a. Therefore, the Dolan-Grady relations arise in this model as well. Again, 
like in model (|2.6|) . for R±± = the equations of motion for the components X± take the 
form of basic commutation relations of the single-mode parafermion. 
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3 Generalized deformed oscillator 



3.1 Main definitions 

Another basic ingredient we shall use in what follows is the generalized deformed oscillator 
(GDO). Therefore, in this subsection we present a brief description of the construction (see 
also Refs. [El ESI). 

The generalized deformed oscillator is defined by the algebra generated by the set of 
operators {I, z, z, N} satisfying the relations 

[N, z] = -6z } [N, z] = 6z } zz = F(N), zz = F(N + 9), (3.1) 

where N is the number operator and the structure function F(.) is an analytic function 
with the properties F(0) = and F(n9) > 0, n = 1, . . .. The structure function F(.) is 
characteristic to the deformation scheme. Here we have introduced a positive real parameter 
9, which will serve later on as a measure of noncommutativity. 
The operators z, z obey the commutation relation 

[z,z]=9A + F(N). (3.2) 

By symbols A± we denote the forward and backward difference derivatives 

A+/(JV) = W + V-fW A _ f(N) = fW-M-V (3 . 3) 

The generalized deformed oscillator algebra can be naturally represented on the Fock 
space of the eigenstates of the number operator N, N\n) = n9\n), (n\m) = 5 nm , n = 1, . . ., 
with the vacuum state defined as z\0) =0. The eigenstates are given by 

1 n 
\n)= r|0), F(n9)\ = l[F(k9). (3.4) 

The operators z and z are the annihilation and creation operators of the oscillator algebra, 



z\n) = y/F(n9)\n-l), z\n) = ^F{{n + l)6)\n + 1). (3.5) 

In what follows we shall denote the Fock space defined by ()3.4j) . ()3.5j) as Hf- On this space 
the operators z and z are mutually adjoint, z = z\ 

Note that when the characteristic function obeys the condition 

F((p+l)0)=O (3.6) 

for some p G N, the creation-annihilation operators satisfy the nilpotency relations z p+1 = 
z p+1 = 0. This means that the corresponding representation is (p + l)-dimensional and we 
have a parafermionic type system of the order p. 
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3.2 The Dolan-Grady relations in terms of GDO 

The simplest realization of the Dolan-Grady relations (j2.1)l can be obtained in terms of the 
generalized deformed oscillator: 

Z = z, Z = z. (3.7) 

As a consequence of the natural grading of the GDO algebra defined by the operator N, for 
this representation of the generating elements Z and Z it is possible to realize the Dolan- 
Grady relations only in its contracted form, u = [9 , 

[Z, [Z, [Z, Z]]]=0, [Z, [Z, [Z, Z]]]=0. (3.8) 

On the other hand, these relations lead to the restriction on the structure function in the 
form of the finite-difference equation: 

A 3 + F(N) = 0. (3.9) 

Since the structure function has to vanish at zero, the general solution of the equation (|3.9|) 
is 

F(N) = Cl N 2 + c 6N, (3.10) 

where ad interim cq,c\ £ R. Three different cases corresponding to the solution (j3.10|) should 
be distinguished: 

1. The case c\ = 0, c ^ corresponds to the usual oscillator given by the Heisenberg 
algebra. The nonlinear supersymmetry for such a system was discussed in Ref. [7]. 

2. For C\ > the polynomial (j3.10|) is a structure function when the parameters obey 
the inequality C\ + Cq > 0. Then the generalized oscillator realizes the half-bounded 
infinite-dimensional discrete series of the unitary representations of the sl(2, R) algebra. 
The generators are, correspondingly, L = 9A + F(N), L_ = z, L + = z, or L = 
—9A + F(N), L_ = z, L + = z for the series V+ or T>~ |Uj with a = (c + c 1 )/{2c 1 ). 

3. In the case c\ < 0, for the solution (j3.10|) to be a structure function one has to impose 
the additional condition ()3.fi|) . This leads to quantization of one of the parameters, 
c o = —ijP + l)ci, and results in the {p + l)-dimensional representation of the su(2) 
algebra with the generators L_ = z, L + = z and L 3 = 6A + F(N). This is the single- 
mode parafermion case. 

4 Noncommutative quantum mechanics and DGR 

Let us consider a 2D quantum mechanical system with noncommutative complex-like coor- 
dinates z, z obeying the GDO algebra (j3.1|) realized on the Hilbert space TCf- The operator- 
valued functions ^(z, z) on Ti F with a finite norm with respect to the scalar product 

($,*) = Tr« F ($t(z,s)*(M)) (4-1) 
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will be treated by us as the states of the noncommutative system. Here the symbols Tr-^ F 
and t denote the usual trace and Hermitian conjugation operations on the space Tip (|H.4j) . 
The operator-valued functions with the finite norm themselves form a Hilbert space, which 
we shall denote as Hp- 

The corresponding derivative operators on the space Tip can be defined as 

0=-0 -1 adz, 9 = -1 adz. (4.2) 

These operators possess the usual conjugation property = —d with respect to the scalar 
product (|4.ip . It is worth noting that the derivative operators do not commute if the structure 
function F(N) is different from a linear one. Indeed, [d, 8] = 6~ l ad A + F(N). 

Let us introduce a U(l) gauge interaction associated with an external magnetic field. 
The noncommutative system admits the following two types of gauge transformations: 

z) = U{z, z)V{z, z), V'\z, z) = tf(z, z)U\z, z), (4.3) 

or 

z) = U(z, z)^(z, z)U\z, z) (4.4) 

with U(z, z) being a unitary operator. The scalar product (j4.1j) is invariant with respect 
to the both transformations. In the commutative limit, 9 — > 0, the transformation (|4.3jl is 
reduced to the usual £7(1) gauge transformation while (|4.4J) becomes trivial. 

We realize the generating elements from (|2.4jl as the gauge covariant derivatives 

Z = d + p[A(z,z)], 2= -B + p[A(z,z)], (4.5) 

where the operator p[0), O G {A(z, z), A(z, z)}, acts on ^ as 

p[0]m = OV, p[0}¥ = 

manifesting the fact that the states \& and carry the charges of the opposite sign. The 
operator-valued functions A and A = A* are the components of the noncommutative vector 
gauge potential corresponding to the case of stationary magnetic field, and we assume that 
the gauge coupling constant ("electric charge") and the mass parameter are equal to unit, 
e = m = 1. Under the both gauge transformations ()4.3j) and f)4.4j) the components of the 
noncommutative vector potential are transformed in the usual way 

A = UAU ] + UdU\ A = UAU ] - UdU ] . 

The commutator of the generating elements has the form 

B = [Z, Z] = - r 1 ad A + F(N) + B(z, z). (4.6) 

We will treat 

B(z,z) = 8A + 3A+ [A, A] (4.7) 

as an external quasi-magnetic field since in the commutative limit the commutator in (|4.7jl 
disappears and it turns into a true magnetic field. Nevertheless, in the noncommutative case 
the quantity ()4.7|) transforms as a connection under the gauge transformations. Actually, if 
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one considers the operator N as the third dependent coordinate, then the quantity (|4.7|) is 
the iV-th component of the gauge connection. 

Let us suppose that the magnetic field depends on iV only. This corresponds to an axially 
symmetric field in the commutative limit. The following choice of the gauge, 

A(z, z) = zf(N), A{z, z) = f(N)z, (4.8) 

where /(.) is a real function, guarantees such a dependence of the magnetic field. Indeed, in 
the gauge ()4.8|) the quasi-magnetic field ()4.7|) acquires the form 

B(N) = A + (F(N)f(N -6) (2- 6f(N -6))). (4.9) 

Now let us analyze the restrictions which follow from the Dolan-Grady relations. In the 
representation ()4.5|) the l.h.s. of Eq. (j2.1j) can be rewritten as follows 

ZAs.(|2ID = - 0' 1 ad(z 2 A^F(iV)) + z 2 (1 - Bf{N + 6)) (1 - 6f{N)) A 2 + B{N) 

+ z 2 (f(N) + f(N + 6)(l- 9f(N)))A 3 + F(N). ^ 

The comparison of the first terms in (|4.fi|) and ()4.1()j) leads to the conclusion that the Dolan- 
Grady relations (|2.1|) with uj ^ can be satisfied only when the structure function of the 
GDO algebra is linear. On the other hand, for u — the expression (|4.10j) vanishes only 
when the structure function obeys the equation (|3.9J) . This leads to the following constraint 
on the quasi-magnetic field: 

A 2 + B(N) = 0. (4.11) 

In other words, the quasi- magnetic field has to be at most linear in N. The function f(N) 
appeared in the gauge ()4.8|) can be found from ()4.9|) taking into account the constraint (|4.11|) . 

The classification of the GDO algebra with the quadratic structure function ()3.10|) was 
presented in the previous section. In accordance with it, the noncommutative coordinates 
z, z and the operator N form the su(2), sl(2,R), or the Heisenberg algebras. Therefore, 
the noncommutative system given by the Hamiltonian (J2.4j) in the representation (J4.5|) can 
be treated, correspondingly, as the system on the fuzzy sphere Sg, on the noncommutative 
hyperbolic plane iff, or on the noncommutative plane IRj. 

Note that since the structure function is quadratic, the simplest solution for f(N) is 

f(N) = f = const. (4.12) 

Indeed, one can verify that in this case the quasi-magnetic field f!4.7|) is a linear function, 

B(N) = f(2 - 6f){2c x N + (c + c x )0), (4.13) 

and hence, it obeys the equation (j4.11j) . Below we shall consider only the solution ([4.12)1 . 
()4.13|) . Later on we shall show that for the given solution the expression (J4.13J) coincides 
with a covariant definition of the magnetic field in the gauge (|4.8|) . 

From the explicit form of the magnetic field (|4.13|) it follows that the parameter / has 
two special values / = and / = 26~ 1 . The first case is rather trivial. It corresponds to the 
model without gauge interaction (a free model (|2.4|) on the corresponding noncommutative 
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space). In the second case the magnetic field ()4.13|) vanishes. Nevertheless, the gauge 
interaction still is present in the model since the corresponding connection (|4.8|) with f ^ 
is not trivial. This resembles the famous Aharonov-Bohm effect but its origin roots here 
in the noncommutativity of the configuration space. As we shall see, there is also another 
critical value of the parameter, / = 0" 1 . For this value the quasi-magnetic field is nontrivial 
while the system is degenerate. Such a case, in turn, resembles the situation with existence 
of the critical value of the constant external magnetic field being coherent with the value of 
the parameter of noncommutativity (see, e.g. |2"T| I22|). 

In the next section we shall discuss the spectral problem for the noncommutative system 
(I2.4J1 on the fuzzy sphere and noncommutative hyperbolic plane. 

5 Reduction to finite-difference QES systems 

Here we show that the noncommutative system given by the Hamiltonian ()2.4j) is quasi- 
exactly solvable for some values of the coupling constant. It is demonstrated by a reduction 
of the Hamiltonian to one- dimensional finite-difference QES operators. 

First of all, one can notice that in the gauge (|4.8jl the operator ad N commutes with 
the Hamiltonian. One can refer this property to the existence in the system of an "axial" 
symmetry. Therefore, the eigenstates of the Hamiltonian can be represented in the form 

V L (z, z) = z m ij(N), or * R (z, z) = ^(N)z m , (5.1) 

where m G Z + . Evidently, these operator functions are eigenstates of the operator ad N. We 
shall call the functions \& l and r the left and right modes. 

Now it is clear enough that the corresponding 2D spectral equation with the eigenstates 
of the form (|5.1|) can be reduced to a one-dimensional problem. But unlike the commuta- 
tive case, here the corresponding ID system is represented by a finite-difference equation. 
One notes that such a situation is typical for any noncommutative system with the "axial" 
symmetry. 

The Hamiltonian (J2.4j) is associated with the nonlinear holomorphic supersymmetry for 
integer values of the coupling constant. The zero modes of the supercharge Q form an 
invariant subspace of the Hamiltonian. We will use their functional form as an anzatz for 
investigating the spectral problem with arbitrary coupling constant. 

Let us temporarily set A = n G N. For the left mode wave functions represented in the 
factorized form 

* L (z,z)=z m V >(N)<l>(N), 

with <f>(N) obeying the equation 

(A + + f)cf>{N) = 0, (5.2) 
the zero mode equation Z n ^/ L = is reduced to 

(l_ W A>(iV)=0. 

For Of 7^ 1 this equation gives the polynomial solution <p(N) = P n _i(N), where Pk{-) denotes 
a polynomial of the k-th degree. For Of = 1 there are no zero modes at all, i.e. there are no 
bounded (normalized) states in the system. 
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The existence of the zero modes means that the Hamiltonian is quasi-exactly solvable 
when the coupling constant A is a nonnegative integer. Now we investigate the question on 
existence of the zero modes for an arbitrary value of the coupling constant. For the left 
modes the Hamiltonian can be reduced, H\^l(z, z) — > H Xm ip(N), to the one- dimensional 
finite-difference operator 



(5.3) 



Km = (*/ - 1) (*W A + A - + (m + 1)HN) A + ) + fF(N) 
+ -/((2m + A + 1)6 f - 2(m + X))A+F(N) + const 

with the function J-{N) defined by 

F(N) = A + F(N) + ^AlF(N). (5.4) 

For Of = 1 the reduced Hamiltonian ()5.3|) is a multiplicative operator. The generating ele- 
ments of the form ()4.5|) are multiplicative operators as well, Z^> ~ ^z, Z^> ~ ^>z. Actually, 
this representation of the Dolan-Grady relations becomes equivalent to the representation 
(13. 7|) . However, the Hamiltonian (|2.4|) acting on the Hilbert space Hp has a trivial dynamics 
in this case. Therefore, later on we shall always assume 9f ^ 1. 

The leading term of the "potential" in ()5.3|) is proportional to the structure function 
F(N). We represent the function ip(N) in the factorized form 

V(iV) = <p{N)<!>(N), (5.5) 

where the function <f>(N) is supposed to obey the equation 

(A + + 7 )0(iV) = (5.6) 

with a real parameter 7. This equation is a slight modification of Eq. (|5.2|) . The factorization 
(J5.5|) allows us to discard the leading term in the equation for <p(N) when the constant 7 
obeys the algebraic equation 

7 2 / 2 



1-07 1-0/' 

that gives a possibility to look for the solutions <p(N) in a polynomial form. 
The first evident solution of the equation (J5.7j) 



(5.7) 



7 = / (5.8) 
corresponds to the zero modes found above. Indeed, the reduction 

HlrMN) - H' x>m <p(N) 
results in the one-dimensional finite-difference operator 

H' x>m = - F(N) A + A_ - (m + 1)(1 - 9f) 2 F(N) A + + /(2 - Of) T x + const. (5.9) 

Here, the function J-(N) is defined by (|5.4j) . while the operator T\ is given by 

T A = F(iV)A+ - AA + F(A^). (5.10) 
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Let Vk(x), k G Z + , denotes the space of all the polynomials in x up to the order k. 
On such spaces the finite-difference derivatives act similarly to the usual derivative, A± : 
Vk(N) — > Vk-i(N). Using this fact and that the function F(N) is quadratic while J-(N) is 
linear, one can notice that the spaces Vk{N) are invariant subspaces of the operator ()5.9|) 
except for the third term. On the other hand, the space V n {N) is an invariant subspace 
of the operator T\ for any A = n 6 Z + . As a result, for nonnegative integer values of the 
coupling constant the operator (|5.9|) is quasi-exactly solvable [T2] . 

The second solution of the equation (|5.7|) is 

As a consequence, the reduced Hamiltonian acting on ^p(N) reads as 

Km = -(! - #/) 2 (*W A + A_ + (m + l)JF(iV) A+)-/(2 - 0/) T_ A _ 2m + const. (5.12) 

Therefore, this operator is quasi-exactly solvable if one imposes the condition A + 2m 6 Z_ 
on the coupling constant. 

Now we pass over to the discussion of the right modes. Unfortunately, in this case we 
failed to find the zero modes of the supercharge Q in a closed form. Therefore, we shall 
investigate the conjectural partial algebraization of this part of the spectrum from viewpoint 
of the reduction of the Hamiltonian. 

For the right modes (j5.1j) the 2D Hamiltonian is reduced to the finite-difference operator 



H* m = (6f - 1)(F(N) A + A_ + (m + 1) F(N) A + ) + f 2 F(N) 
+ -/((A + 1)9 f + 2(m - A)) A + F(N) + const. 



(5.13) 



This Hamiltonian can be obtained from ()5.3|) by means of the formal change A — > A — 2m. 
Using this property, one can reproduce all the corresponding results for the right modes. 
The factorization ()5.5j) with the conditions (|5.fij) and ()5.7j) for the solution ()5.8j) leads to the 
finite-difference operator containing the term proportional to the operator T\_ 2m - Hence, 
for A — 2m £ Z + the polynomial space V\-2 m {N) is an invariant subspace of the obtained 
operator, i.e. it is QES operator. For the second solution (|5.1H) the reduced Hamiltonian 
contains 71 A , and so, it is QES operator for non-positive integer values of the coupling 
constant. 

All the reduced finite-difference operators, (|5.9jh (|5.12|) and those for the right modes, 
can be represented as QES operators of the form (jA.lOj) with the following coefficients: 

1. The solution (Oil : 

ai = 0, a 2 = 0, bX= Of) 2 , 6+ = c , b± = 0, 6; = |, 

foo = - (c + 2(m + l) Cl ) (1 - 0/) 2 , 61 = - 0(m + 1) (c + (m + l) Cl ) (1 - Of) 2 , 

where N = AorN = A — 2m for the left or right modes, respectively; 
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2. The solution fl5Tj| : 



ai = 0, a 2 = 0, 6+ = - 6+ = c (1 - Off, b± = 0, 

b- = j(l- 9f) 2 , &o = - c - 2(m + l) Cl , 61 = - 0(m + 1) (c + (m + l) Cl ) , 

where N = — A — 2m or N = —A for the left or right modes, respectively. 

It is worth noting that for Of = 2 the both operators ()5.9)1 and (j5.12j) are exactly solvable 
since the term proportional to T\ vanishes. 

Let us briefly discuss a normalizability condition for the found solutions. Since the 
function <p(N) is a polynomial, the normalizability has to be provided by the factor 4>(N). 
The solution of the equation ()5.6|) has the form <p(N) ~ (1 — 0j) N . Therefore, the norm of 
the corresponding state \& £ H F can be represented as 

oo 

11*11*, = Tr WF (* f *) = Cvj>*(l - fl 7 ) 2n + • • ■ ■ (5-14) 

n=l 

Here the dots denote less divergent terms in the sense that if the first term in (|5.14jl is conver- 
gent, then the others are convergent as well. The factor CV is related to the normalization 
constant and the integer parameter K is defined by the type of solution. The requirement 
of convergence of the series leads to the condition 

|1-07|<1. (5.15) 

Therefore, we have to impose the following restrictions on the values of the parameter /: 
Of £ (0, 2) for the solution (pTSjl . and Of £ (-oo, 0) U (2, oo) for the solution (l5~TTJ) . 

Of course, in the case of the fuzzy sphere the norm (|5.14|) is finite because of the finitness 
of the corresponding Hilbert space, and no additional condition appears for it. 



6 Intrinsic algebra and quasi exact solvability 

In this section we show that the noncommutative system under consideration really is quasi- 
exactly solvable for any value of the coupling constant. The demonstration is realized in 
terms of an intrinsic algebra. 

6.1 Intrinsic algebra 

In general, the operators Z = Z and Z = Z together with the contracted Dolan-Grady 
relations ()3.8|) recursively generate the infinite-dimensional contracted Onsager algebra [§]: 

[Zk, %i\ = B k+ i + i, [Zk, B[] = Z k+i , \B k , Zi\ = Z k +i, 

[Z k , Zi\ = 0, [Z k , Z t ] = 0, [B k , B t ] = 0, 

where k, I £ Z + and Bq = is implied. This algebra can be extended discarding the last 
condition, i. e. by treating Bo as a nontrivial generator. Then the extended algebra will be 
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generated by the set Z, Z and B . The element B splits up the contracted Onsager algebra 
into the three subalgebras of the grade +1, and —1, given, correspondingly, by the span 
of Z m , B m and Z m . In other words, the generator B is a grading element which provides 
the extended algebra with the triangular decomposition. This grading is different from that 
considered in Ref. 0. 

The algebra ()6.1j) admits the infinite set of the commuting quadratic charges 



k « k-1 



J a = g £ { Z p~ii Z k- P ) ~ g B p B k-P ~ (6-2) 
p=l p=l 

with G Z + , which contains the Hamiltonian (|2.4jl . = 2_ff A , and the grading operator, 
J° A ~5o. 

In the gauge ()4.8|) the generating elements of the form ()4.5|) read as 

Z= -^adz + Zpff], Z = -r 1 adz + /p[^]. 

These operators generate the following algebra intrinsic to the noncommutative system: 

[Z, Z] = 2 Cl (S + f(2-8f)G), 
[Z, P] — 2ci (1 — Of) G, 
[P, Z] = 2 Cl (l- 9 f)G, 

[Z,G] = (l-9f)V, [S. V] ■ P. (0.3) 

[G, Z] = (l-0f)P, 
[Z, 5] = Z, 

[S, G] = [Z, P] = [Z, P] = 0, 

where the operators P, P, G and S in the chosen gauge are 

V = z, V = z, G = N + 9 Co + Cl , S=-6~ 1 adN. (6.4) 

2ci 

The algebra has the following two Casimir operators: 
d = {P, P} - 2 Cl G 2 , 

G 2 = {#Z + (1 - 9f)V, 9Z+(1- 9f)V) - 2ci (05 + G) 2 



[P, P] 


= -2 Cl #G, 


[G, P] 


= 0P, 


[P, G] 


= &D, 


[S, P] 


= v, 


[P, 5] 


= v, 


[5, Z] 


= z, 



(6.5) 



The algebra (Q is sI(2,R) ©sI(2,R) ( Cl > 0) or su(2) © su(2) (ci < 0). In the chosen 
gauge it can be shown by representing the generators of the algebra in the form L z , L 2 , 
and R z , R z , Rn, where L a and R a denote the left and right multiplications. In the given 
representation the Casimir operators ()6.5|) are equal to the same number, 

Ci = C 2 = ^ (c\ - cl) . (6.6) 

It is necessary to note that the values 9f = 0, 1, 2 are special from the viewpoint of 
the algebra ()6.3|) . Although the algebraic content of it remains the same, the algebra is not 
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generated by the operators Z and Z, i. e. it is not intrinsic any more. In these cases the 
generating elements produce the sl(2,R) or su(2) only. 

The operators £>, T> and G can be treated as dependent covariant coordinates of a non- 
commutative membrane. Indeed, the equivalent set of coordinates 

V + V V — V 

Xi = — / , n X-2 = i — , , , , X% = G 

x/sR \/8N 

obeys the condition 

X 2 + X 2 - eXl = e9 2 a (1 - a) = -er 2 , (6.7) 

where e = sign ci and a — {c\ + co)/(2ci) (ci > 0) or a = —p/2 (ci < 0), p G Z+. So, 
the operators Xj form the set of noncommutative covariant coordinates of the fuzzy sphere 
Sq (ci < 0), or of the noncommutative hyperbolic plane Hg (ci > 0). The operators Z, 
Z and 5*, in their turn, can be interpreted as the covariant derivatives. Note that in the 
noncommutative system under consideration, unlike the constant magnetic field case [22J, 
the covariant coordinates are not proportional to the covariant derivatives. 

The intrinsic algebra gives us the possibility to define the magnetic field in a covariant 
way. In general, if the derivatives form a Lie algebra, [<9j, dj] = fij k dk, then the corresponding 
gauge covariant field strength can be defined as iF^ = [D^, Vj] — fij k V k (e. g., see Ref. [22j), 
where T>i are gauge covariant derivatives. From this definition and the intrinsic algebra it 
follows that in the present case only the "transversal" component (i = 3) of the magnetic 
field, Bi = ^f^ k Fj k , does not vanish. In the gauge ()4.8|) this component exactly coincides 
with the definition (J4.13)) . It is interesting to note that like in the commutative case 
the magnetic field is a linear function of the "transversal" coordinate. 

In terms of the generators of the algebra ()6.3|) the higher operators of the contracted 
Onsager algebra are 

z m = (2 Cl ) m (z + e-\i - - (1 - 6f) 2m )v), 

B m = (2 Cl ) m (S + 6-\l - (1 - 6f) 2m )G), 

where 9 f 7^ 0, 1, 2. Applying this form of the generators of the Onsager algebra to the infinite 
set of the commuting charges ([6.2)1 . one can obtain 

Ja = 0^2 -Of) ( 2 (1 " (1 " 6f)2k) Hx ~ Cl ((1 " 61)2 " (1 " ° f)2k) S{S ~ A) ) + ■ ■ ■ ' 

where dots denote a linear combination of the Casimir operators (|6.5jl . This representation 
leads to the conclusion that the Hamiltonian (|2.4jl and the operator S form a complete set 
of nontrivial integrals of the system under consideration. 



6.2 The commutative limit 

Our aim is to apply the found intrinsic algebra of the noncommutative system with the 
Hamiltonian ( I2.4p for investigation of its spectrum. The algebraic scheme we shall discuss is 
an analogue of the Fock space construction. It is based on the introduction of the "ground" 
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vectors in the Hilbert space of the system (|2.4jl . These vectors have special properties with 
respect to the action of the generators of the intrinsic algebra. 

First of all, let us note that in the commutative limit 9 — > the algebra (|6.3|) is converted 
into the intrinsic algebra of the system (|2.4|) living on the (pseudo) sphere [TI]. Indeed, with 
the rescaling 

G -£~r c '- 2ft 

where the parameter (3 is related to the radius of the (pseudo)sphere, the algebra ()6.3j) 
changes into the intrinsic algebra appeared in Ref. ^T] in the limit 9 — ► 0. 

Since the commutative case is much simpler than the noncommutative one, we shall start 
the discussion of the covariant algebraic construction from it. 

We introduce a vector ( "ground" state) of the Hilbert space which obeys the conditions 

Ztf<°> = Va(G)^°\ ZV<® = V(3(G)¥°\ S^ 0) = s¥°\ (6.8) 

Here a(.) and j3(.) are some functions, which will be defined later, and s G R. We start 
with the real eigenvalues of the operator S since no quantization condition follows from the 
intrinsic algebra. Besides, we know ^T] that the eigenvalues do are real in general, since the 
magnetic flux gives a contribution to them and for a noncompact surface it is not quantized. 

From the commutation relations for the generating elements Z and Z the following con- 
straint on the functions a(G) and [3(G) appears: 

1 G 2 -) gi(G) + Gq{G) = G- 4/5 2 s, (6.9) 



2 

where g(G) = a{G) + (3(G). The solution of the differential equation is 

0(G) = 1 - ^ - ^-j. (6.10) 

The parameter m is the integration constant. One should regard this constant as an integer 
since it corresponds to the quantum number of the planar angular momentum of the system 

EH- 

We will look for the eigenstates of the Hamiltonian in the form 

$ = a(G)¥ 0) , (6.11) 

where we suppose that cr(G) is a polynomial. Then the Hamiltonian ()2.4j) is reduced to the 
second order differential operator 



H' x =- 4,f- (C 2 - |) + ((g 2 - |) (a(G) - 3(G)) - SB'g) 
+ G( 2/ + ,- 1)+ (G-f)( a ,(G, + «) 



d 

dG 



(6.12) 

const 



16 



acting on the function cr(G). 

Let us fix the residuary arbitrariness in the definition of the functions a(G) and (3(G) re- 
quiring that the operator (|6.12|) is quasi-exactly solvable. Therefore, it should be represented 
in terms of the generators of the st(2, K) algebra 



J+ = 0>± _ (N _ 1)G , J, = G±-^1, J_ = ±, (6,3) 

where if the real parameter N is a nonnegative integer, it defines dimensionality of the 
corresponding finite-dimensional representation. Moreover, the operator (J6.12)) is a quadratic 
polynomial in the generators (|6.13|) since it is a second order differential operator [T2j . 
Eventually one can conclude that the operator (j6.12|) can be represented in the form 

T = c 2 /3 2 J 2 - A(3 2 Jq + a + J + + a J + a_ J_ + const, (6.14) 

where the parameters a±, ao and N have the following eight solutions: 

|ci (ci + 16/3 2 (m - s)) , 
| Cl ( Cl - 16/? 2 s), 
|ci(ci - I6f3 2 (m - s)), 
| Cl ( Cl + 16/3 2 s), (6.15) 

where e = ±1. So, if the parameters of the system take such values that N £ Z + for any of 
the solutions, then the system is quasi-exactly solvable. Any QES part of the spectrum is 
a superstructure over the corresponding "ground" state. Therefore, one can interpret these 
states as some QES "vacua" of the system. Here and in what follows we will understand 
under the QES "vacuum" a state over which one can construct a certain number of states 
of the QES part of the spectrum following an oscillator-like procedure based on an intrinsic 
algebra. 

The number of the sl(2, M) representations ()6.14|) coincides with the number of the solu- 
tions found in Ref. [TT]. At the same time, it is worth noting that the method based on the 
intrinsic algebra is completely covariant since all the generators of the algebra are covariant 
quantities. 

We would like also to remind that this commutative system possesses a discrete symmetry 
Z 2 ® Z 2 ® Z 2 , which intertwines the solutions (|6.15|) between themselves. The detailed 
discussion of the symmetry can be found in Ref. |llj . 



N 


= eA, 




a + 


= e, 


a 


= 4e/3 2 (2s - A) , 


a_ = 


N 


= e(A- 


m), 


a + 


= e, 


«o 


= 4e(3 2 (2s-m- A), 


a_ = 


N 


= e(A- 


2s), 


a + 


= e, 


«o 


= -4eA/3 2 , 


a_ = 


N 


= e(A + 


m — 2s), 


a + 


= e, 


a 


= 4e(m - A)/? 2 , 


a_ = 



6.3 The noncommutative case 

To investigate the spectral problem for the noncommutative system with the Hamiltonian 
(12. 4|) in terms of the intrinsic algebra ()6.3|) . we introduce a "ground" state obeying the 
conditions ()6.8|) but with the only change 

S¥ 0) = m¥°\ meZ. (6.16) 
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One can treat this number as an "azimuthal" quantum number since it corresponds to 
the existence of "axial" symmetry in the system with operator S as its generator. We 
assume this in accordance with the representation (J5.1|) of the eigenstates. Therefore, in the 
noncommutative case the conserved quantity S is an analogue of the angular momentum 
operator. We would like to note that although we use formally the same conditions ([6.8)1 . 
one should remember that the operators T>, T> and G do not commute and the ordering of 
the operators is important. 

From the commutation relations of the generating elements ()6.3j) and the definition of 
the "ground" state (j6.8|) it follows that the functions a(G) and (3(G) are not independent 
but obey the constraint 

2c x G((l - Of) {a{G) + 13(G)) + 6a(G)(3(G)) + VV(l-6f + 9a(G)) A + (3(G) 

+ VV(l-6f + 6f3(G))A_a(G) = 2c x (m + f (2 - Of) G) . (6.17) 

This is the noncommutative analogue of the restriction (|6.9|) . 

For investigation of the spectral problem it is convenient to represent the Hamiltonian as 
the linear in Z, Z form: 

H x = - 26- 1 (1 - Of) (VZ + VZ) + 2c 1 r 1 / (2 - Of) G (G - X9) 

+ a Ci 9- 1 gs + 2 Cl s (s - A) - e- 2 d (i - eff + e- 2 c 2 . 

Using the definition of the "ground" state ()6.8|) . the spectral problem for the Hamiltonian 
in this form is reduced to that of the finite-difference operator 

H A = 2 [Of - 1) (vf>(l - 6f + 6(3{G)) A+A_ + (VVa{G) - VV(3(G) 

+ 2 Cl (1 - 6f) G)A + + r 1 (VVa(G) + VV(3(G)j 
+ 2cJ0- 1 (2 - Of) G(G - X6) + AdmO^G + const (6.19) 

acting on the function u{G). Here the operators A± are defined as the difference derivatives 
(13. 3|) but with N changed for G. The operator (|6.19|) is a finite-difference operator in the 
variable G since T>T> and T>T> are the quadratic polynomials in G, 



2 J 4ci V 2 J 4ci 

as this follows from the intrinsic algebra ()6.3)1 and the form of the Casimir operators ()6.5)1 , 
()6.6|) . The evident advantage of the representation ()6.18|) is that the finite-difference operator 
()6.19|) is linear in the functions a{G) and (3(G). 

Representing the functions a(G) and (3(G) in the form 

a(G) = m, m = sb, 

vv vv 

one can see that if the functions a(G) and (3(G) are polynomials, then the operator ()6.19|) 
is a finite-difference operator with polynomial coefficient functions and vice versa. 
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Let us fix the form of the functions a(G) and j3{G) by demanding the Hamiltonian (|6.18|) 
to be a QES operator of the form (jA.10|) for some values of the parameters. Moreover, we 
require that the functions to obey the constraint ()6.18|) . In comparison with the commuta- 
tive case the noncommutative one is drastically different. There are four solutions for the 
functions a(G) and (3{G) providing the quasi exact solvability of the Hamiltonian (j6.18J) . 
The explicit form of the function is not important. The corresponding solutions for the 
coefficients in the representation (jA.lOJ) are 

ai,2 = 0, b%= T ^1(1-0/)^ 6± = 2c 1 (2j+/i((l±l)m-A))(l-^) 1 ^\ 

bt=± e{Cl ~ Ci eC0) («* - (4j + 2/i((l ± n)m - A) + l) Cl ) (1 - Of)***, (6.20) 

where e = ±1 and fi = ±1 parametrize the set of solutions. The main difference of these 
solutions from those of the commutative case f)6.15j) is the absence of restrictions on the 
parameter N = 2j + 1. Therefore, the operator (j6.19|) is quasi-exactly solvable since we can 
always choose the parameter N to be natural. 

The functions a(G) and (3{G) corresponding to the solution (|6.2U|) depend on the param- 
eter j. Therefore, one can consider that it labels the "ground" states as well, 

Here we recovered the dependence of the state on the quantum number m (16.16)) . Hence, 
unlike the commutative case, the noncommutative system admits the infinite number of 
the "ground" states for every value of the quantum number m. These states serve as QES 
"vacua" of the model. 

Like the commutative case, the noncommutative system has a discrete symmetry. Indeed, 
evidently that the substitution 

Co -> -c (6.21) 

converts the solutions ()6.20|) with e = 1 to those with e = — 1 and vice versa. This transfor- 
mation is similar to the transformations of the commutative case since it converts one QES 
operator ()6.19|) into another of the same type. Besides, in this case there is another discrete 
transformation 

9 -»• -9, /-►-/, A -> -A, m -> -m. (6.22) 

This mapping induces the following change of the parameters (|6.2(J|) and the operators (|A.5|) : 

where a — 0, ±. The transformation of the generators is provoked by the change of the 
sign of the parameter 9. This follows from the definitions of the forward and backward 
derivatives ()3.3|) and of the generators (jA.5|) themselves. In this way, the transformation 
(I6.22J) intertwines the QES operators corresponding to the solutions ()6.20|) with /x = 1 and 
/i = — 1, respectively. The whole group of the discrete symmetry consists of four elements 
including the identity mapping. The transformations (|6.21|) and (j6.22J) are associated with 
the generators of a finite generated Abelian group (22] • This discrete group is Z 2 (g> Z 2 since 
all its elements are involutive. So, like in the commutative case, this discrete symmetry is a 
kind of duality. 
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7 Discussion and outlook 



Let us summarize briefly the obtained results and discuss some open problems that deserve 
further attention. 

The Dolan-Grady relations provide a rich algebraic structure, which arises in numerous 
branches of mathematical physics. These relations originally appeared in the framework 
of various integrable models [TJ |21 El HI E] Later it was observed that the Dolan-Grady 
relations are necessary conditions of the anomaly- free quantization jHJ of a pseudo-classical 
ID systems with the nonlinear holomorphic supersymmetry [7|. This observation revealed 
the intimate connection of the relations with the QES systems. The algebraic origin of the 
construction allowed us to apply it to the two-dimensional systems [TUl ITT] . In the present 
paper we demonstrated that 

• The Dolan-Grady relations arise in low- dimensional matrix models. 

The matrix models naturally lead to appearance of the noncommutative geometry [r51IT4"ll2*4"] . 
Using this observation as a hint, we applied the algebraic scheme described in Ref. P to 
the two-dimensional systems on a noncommutative background subjected to an external 
U(l) gauge interaction. We started by postulating that the coordinates form the generalized 
deformed oscillator algebra jTHj . The Hamiltonian of the system was constructed in terms of 
generating elements realized as gauge covariant derivatives on the noncommutative manifold. 
We showed that the generating elements obey the Dolan-Grady relations only when 

• The noncommutative coordinates generate su(2), sl(2, M) or the Heisenberg algebras. 

This allowed us to identify the resulting noncommutative manifold, correspondingly, with 
the fuzzy sphere S%, the noncommutative hyperbolic plane Hq or the noncommutative plane 
M.g. The cases of su(2) and s[(2, M) were discussed in detail. 

We showed that the parameter / defining the intensity of the quasi-magnetic field has 
the three critical values: / = 0, / = 6' 1 and / = 26~ 1 . The first case is trivial since it 
corresponds to the absence of the gauge interaction (B — 0). The second case is analogous 
to the situation for the particle on the fuzzy sphere and the noncommutative hyperbolic 
plane in the presence of the constant magnetic field B = 9~ l [22] (see also Ref. j2JJ for the 
case of the noncommutative plane). For this special value the systems are degenerate. In 
the third case, / = 2(9 _1 , the quasi-magnetic field vanishes while the gauge interaction still 
is present in the model since the corresponding connection is not trivial, / 7^ (see Eqs. 
(J4.8|) . (|4.12j) ). This resembles the famous Aharonov-Bohm effect, but here its origin roots in 
the noncommutativity of the configuration space. Note that there is no analogue of such a 
phenomenon in the constant magnetic field case. 

By the reduction of the Hamiltonian to one-dimensional finite-difference operators, it was 
showed that the constructed noncommutative systems are quasi-exactly solvable for some 
values of the coupling constant. In the commutative limit the systems on the fuzzy sphere 
and noncommutative hyperbolic plane lead to QES systems on the corresponding Riemann 
surfaces [TTj . 

The algebraic content of the given model was revealed by the construction of the intrinsic 
algebra. Using the notion of the intrinsic algebra, in Ref. JTJ we proposed a general algebraic 
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method of constructing the covariant form of the integrals of motion. Developing this idea, 
here we demonstrated how 

• The spectral problem can be treated in terms of the intrinsic algebra. 

The advantage of this approach is its complete covariance. For example, one does not need 
to impose any gauge conditions since the intrinsic algebra operates with covariant quantities. 
The application of this method allowed us to show that the noncommutative systems under 
consideration are quasi-exactly solvable for any value of the coupling constant. 

In the noncommutative case the reduction of the Hamiltonian and the application of in- 
trinsic algebra to the spectral problem lead to some finite-difference equations. To investigate 
their quasi exact solvability, 

• A deformation of the sl(2, R) scheme for finite-difference QES systems was proposed. 

In this approach the sl(2, R) algebra is changed for its nonlinear analogue (jA.6|) . We showed 
that the deformed scheme is not equivalent to the approach based on the usual sl(2, R) 
algebra [23] if one requires for finite-difference QES operators to be the polynomials in the 
generators of the corresponding algebras. 

Surprisingly enough, but the system ()2.4j) in a non-constant magnetic field on the non- 
commutative plane (the coordinates form the Heisenberg algebra) is more complicated than 
the systems on the fuzzy sphere and the noncommutative hyperbolic plane investigated in 
this paper. Therefore, it is of interest to investigate this system and verify if it inherits such 
a property as the quasi exact solvability. Furthermore, it would be interesting to apply the 
observed analogy between the matrix equations and the parafermionic trilinear commutation 
relations [TH] to construct new solutions. Such a consideration will be presented elsewhere 

Acknowledgements 

We thank W. Orrick and J. Perk for useful comments. M.P. is also grateful to A. Tur- 
biner for valuable discussions and kind hospitality during his visit to Universidad Nacional 
Autonoma de Mexico, where part of this work was done. The work was supported by the 
grants 1010073 and 3000006 from FONDECYT (Chile) and by DICYT (USACH). 

A Deformed 5 [(2, R) scheme 

Here we discuss a deformation of the sl(2, R) scheme related to one-dimensional finite- 
difference operators, and compare the approaches based on the deformed algebra and on 
the usual one. 

In Ref. Smirnov and Turbiner proposed the approach to finite-difference QES equa- 
tions based on the following representation of the Heisenberg algebra: 

a = A + , b = x(l-6A_), (A.l) 

where x is a real variable and A± are the forward and backward difference operators (J3.3|) 
with respect to this variable. These operators obey the standard commutation relation, 
[a,6] = l. 
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Like in the standard approach to ID differential QES equations, an operator for the 
corresponding spectral problem is supposed to be a polynomial in the sf(2,R) generators 

J+ = b 2 a - (N - 1)6, J = ba-^-^, J_ = a. (A.2) 

When N G N, the space of polynomials V^-i(x) is invariant with respect to the action of 
these operators. Hence, any polynomial in the generators (jA.2|) is an invariant operator on 
this space as well, i. e. it is a QES operator. 

We are interested in the spectral problem for a three-point finite-difference equation, 

A(x)^(x + 9) + B(x)V(x) + C(x)^(x -9) = E^(x). (A.3) 

Ref. claims that all such QES equations with finite number of polynomial eigenfunctions 
are classified via the cubic polynomial element of the universal enveloping algebra of sl(2, R) 
taken in the representation (jA.l|) : 

T = A + ( J+ + QJl) + A x J 2 (#J_ + 1) + A 2 J J- + A 3 J + A 4 J- + const, (A.4) 

where the coefficients are arbitrary real parameters. Below we shall show that this operator 
is not the most general among QES operators corresponding to the three-point problem 
(|OJ). 

Let us introduce the finite-difference operators which are linear in the finite-difference 
derivatives, 

J$ = x 2 A± - 2jx, = xA T - j, Jt = A=p (A.5) 

The half-integer parameter j (2j G Z + ) is an analogue of "spin" in the usual sl(2, R) scheme 
[T2*] . These linearly independent operators are finite-difference analogues of the st(2,R) 
generators, which can be reproduced in the commutative limit, 9 — > 0. The operators (jA.5|) 
have another useful property of the s[(2,R) generators. They have the invariant finite- 
dimensional subspace, which is the space of polynomials V^jix), 2j G Z + , d\mV2j(x) = 
2j + 1 = N. This allows us to use the generators (jA.5|) for the construction of finite- 
difference QES operators in the manner analogous to the continuous case. The exhaustive 
investigation of this construction lies out of the scope of the paper and here we will discuss 
just some aspects of it. 

The both sets J+ and J~ , fi = 0, ±, form the nonlinear subalgebras of the form 

[J±, Jf] =2J ± ±fl(2j-l)J±, 

[J ± , J}] =J^T0 (2JJ J ± + (2j + 1)(J ± + j)) , (A.6) 
[J±, 4} = J^T T , 

where the operators 

T± = 1 ± 9JT (A.7) 

are the discrete forward and backward translation generators, T±ip(x) = ip{x ± 9). The 
intertwining commutation relations between the sets Jt and J~ are 

[JT , Jj] = (2J T ± 9JT) T±, [J+, J-] = 9 (2J+ J + - J+ + (2j + 1) J") , 

[J3=, J±] = (J± ± C( J T + j)) T ± , [J + , j -] = J + - J ~, 

[Jf, J ± ]=J^, [J+, Jl]=0. (A.8) 
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It is obvious that in the commutative limit, 9 — > 0, the nonlinear commutation relations 
(|A.6J) . (jA.8|) tend to the usual s[(2,R) algebra. Furthermore, for the realization (|A.5|) of the 
nonlinear algebra the following quadratic identities hold: 

44 -j%Jt- 4 ± o (4 + j) J- - j{j + 1) = o, 
o {4 + j) {Jo +j-i) + J+-J$ = o, 
44 - j-j - u + 1) (j+ - j-) = o, 

J-T+ -4 + 2j 9 (J- + j) = 0, (A.9) 
(J 0+3 )T_-4-j = 0, 
J~T_ - J+ = 0. 



Taking into account these identities and the commutation relations (jA.6|) . (jA.8|) . the most 
general three-point operator ()A.3|) . quadratic in the generators (jA.5|) . can be written as 

Hqes = a x 4J- + a 2 J-J +b+4 + b+4 + btJ+ ^ AiQ ^ 
+ b+ + 6g Jq + bZJZ + const. 

By the construction this operator is quasi-exactly solvable. 

The generators of the s[(2,M) algebra in the representation (|A.2|) can be represented as 
polynomials in the operators (|A.5J) . 

J + = (J+-9(4+j))T_, Jo = Jj, J. = JZ. (A.ll) 

Hence, the QES scheme based on the st(2,R) generators (|A.2|) is included into that based 
on the operators (|A.5|) . On the other hand, the operators (jA.5|) can also be represented in 
terms of (jA.2|) via the two last equalities in (jA.ll|) and 



4 = J + + 9 (J + j) {9 (J +j) J- + 2J - 1), Jl = J-(l + 9J-)- 1 , 
J- = J+T+ + 9( Jo + j) , Jo" = Jo + ( Jo + j) J- 



(A.12) 

According to these relations the shift operator T + (jA.7|) is linear in the generator J_. As 
we see, the correspondence is not of a polynomial form for the operator Jj\ Therefore, the 
sl(2, R) scheme (IA.2J) is not equivalent to that based on the operators (jA.5|) if one restricts 
QES operators to a polynomial in the generators of the (|A.2|) and (|A.5|) forms. One can 
also notice that the 5-parametric (up to overall shift) QES operator (jA.4|) is not the most 
general QES operator in the framework of sl(2,R) scheme. Indeed, for the case 63 = the 
QES operator (jA.lOj) can be represented as a 5-th degree polynomial in the generators (|A.2|) 
with 7 arbitrary real parameters. 

In addition, we note that if one imposes the following restrictions on the parameters of 
the QES operator (lA~Tol) 

ai =0, a 2 = 0, b- = -6+ (A.13) 

then it becomes exactly solvable. 
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